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Рассмотрен модифицированный метод 5Т5-САРСН( 1, 1). Модификация заключалась в отказе от предположения о нормальном за- 
коне распределения логарифмов дневных приращений временного ряда и в использовании для их описания ЕтооіЫу Тгипсаіесі 
а-51аЫе (575 )- распределения (гладко усеченного а-устойчивого). Параметры метода найдены методом максимального правдо- 
подобия. Проведено статистическое исследование надежности предложенного алгоритма и показано уменьшение автокорреля- 
ции в структуре данных, использованных для анализа. Метод применен для прогнозирования цен акций РАО ЕЭС РТС с лагом 5. 


1. Введение 

Изучение свойств, вычисление параметров и 
определение вида распределения некоторого сто- 
хастического процесса, лежащего в основе рыноч- 
ных флуктуаций, является центральной задачей 
эконометрики. Знание распределения необходимо 
при конструировании эконометрических методов 
(АКСН, О АКСН, ЕСАКСН, ЕІСАКСН, ЕІЕ- 
ОАК.СН и др., подробнее о методах см. [1]), при 
оценке предельной величины риска ѴАК, при рас- 
четах вероятных в будущем значений временных 
рядов и при определении асимптотического пове- 
дения плотностей функций распределений. По- 
следнее особенно важно, так как редкие события, 
определяющие форму и вид их хвостов, соответ- 
ствуют получению наибольшей возможной прибы- 
ли или несению наибольшего вероятного убытка. 

В преобладающем большинстве случаев лога- 
рифмы дневных приращений котировок финансо- 
вых инструментов (акций, облигаций, свопов, оп- 
ционов и т. п.) не имеют нормального распределе- 
ния [2-4]. Это связано с тем, что у эмпирической 
функции плотности распределения, построенной 
на основе таких логарифмов, существует ненуле- 
вой эксцесс и асимметрия, присутствует вытяну- 
тость функции плотности в е-окрестности точки 
математического ожидания, а также наблюдаются 
так называемые «толстые хвосты», когда вероят- 
ность значительных изменений цен выше, чем для 
нормального распределения. Все эти факторы 
усложняют или делают невозможным применение 
известных эконометрических методов: АКСН(р), 
ОАКСН(р,д) и др., которые изначально были по- 
строены на допущении о нормальном распределе- 
нии приращений и остатков. 

Неудовлетворенность практических участников 
финансового рынка результатами, полученными 
на основе нормального приближения, заставила 
исследователей искать новые распределения и раз- 
рабатывать новые подходы для обработки эмпири- 
ческих финансовых данных. Так, в работах [5-7] 
для описания временных рядов было использовано 
обобщенное распределение Парето, в [8, 9] - обоб- 
щенное /-распределение Стьюдента, в [3] - распре- 
деление Лапласа, в [10] - а-устойчивое распреде- 


ление. Однако в настоящее время все большее ра- 
звитие получает идея комбинирования вышепере- 
численных распределений с нормальным (см., на- 
пример, [11]). Идея заключается в отсечении хво- 
стов исходной функции плотности и в замене их на 
хвосты нормального распределения. 

В данной работе рассматривается модифициро- 
ванная модель ОАКСЩ1 , 1 ). Модификация заклю- 
чалась в отказе от предположения о нормальном 
законе распределения логарифмов дневных прира- 
щений временного ряда и в использовании для их 
описания БтооіЫу Тгипсаіесі а-8іаЫе (8Т8) -ра- 
спределения [11]. Построение ЗТЗ - рас п редел е н ия 
осуществляется нахождением параметров нор- 
мальных распределений, формирующих хвосты, и 
вычислением первого и второго начальных момен- 
тов. Эффективность предложенной модификации 
5Т$-ОАКСН(1,1) показана с помощью имита- 
ционного моделирования цен акций РАО ЕЭС 
РТС. При этом было использовано 371 значение 
долларовых котировок за период с 4 января 2003 г. 
по 30 июня 2004 г. (данные предоставлены компа- 
нией РТС, 1йІр://\ѵ\у\ѵ.гІ8.т). 

2. Общие положения 

Рассмотрим классический метод ОАКСН(1,1) 
[12]. Пусть И„, п= 1,2,..., - некоторый временной 
ряд. Допустим, что справедлива авторегрессионая 
зависимость: 

а; = уѴ + а г/;_, + /Гт;-, = со + аи п _, + /Гст„ 2 _,, ( 1 ) 

где у>0, а>0, /3>0 - некоторые коэффициенты мо- 
дели, К>0 - долговременное среднее отклонение в 
структуре данных, у+а+р = 1, со=уѴ, 
и=\п(Н^—\п(Ь„_^) — логарифмы приращений значе- 
ний временного ряда И„, сг„ - дневная волатиль- 
ность, я=1,2,... 

Вместо ц„=1п(/г„)-1п(/г„_ 1 ) можно использовать, 
например, оценку м„=1п(/г„ — В обоих слу- 
чаях и„ будут зависимыми случайными величина- 
ми. Поэтому для и„ потребуем существования по 
крайней мере двух первых условных начальных и 
центральных моментов. Для корректного вычисле- 
ния последних предположим, что на вероятност- 
ном пространстве (П,Г,Р), где Р - ст-алгебра под- 
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множеств О, задана фильтрация Р=(Р„)„ Л , состоя- 
щая из ст-подалгебр Р„ таких, что Р т аР п сіР, если 
т<п. При этом события из Р„ будем интерпретиро- 
вать как доступную на момент времени (п- 1) ин- 
формацию. 

Пусть далее Е(и,\Р н1 )=а, и РКи,\Р^)=^=Р(^\Р н1 )-а1, 
і= 1,2,... Если предположить, что и~Ща„,сг„), то в 
соответствии с методологией САКСН(1.І ) выпол- 
нено соотношение: 

и„ = ст„<е„ +а п ,п =0,1,2,..., (2) 

где е„~/Ѵ(0,1) - стандартная нормальная случайная 
величина. Это позволяет производить имитацион- 
ное моделирование будущих значений временного 
ряда И„ по вычисленным согласно (1) волатильно- 
стям сг„ + і: 

ЧК + і) = НК) + ѵ„ +і е п+ і +Ч, + Р П = 0,1,2,... 


Р\ = \ р 2 = 2 (1)СІІ. 

—00 Ь 

Для определения р { и р 2 необходимо пользовать- 
ся квадратурной формулой и вычислять соответ- 

а со 

ственно интегралы | % в (х)сІх, | % 0 (х)с1х. 

-со Ь 

Вероятности р^ и р 2 играют важную роль при де- 
терминации параметров нормальных распределе- 
ний с плотностями й,(х), И 2 (х). Справедлива сле- 
дующая теорема. 

Теорема 1. Пусть 8(х,а,Ъ,а,(5,с,р) - функция 
5Т8 -распределения с плотностью Дх), удовлетво- 
ряющей (4). Тогда параметры нормальных распре- 
делений с плотностями А|(х), И 2 (х) вычисляются по 
формулам: 


Однако в произвольном случае, если закон ра- 
спределения и„ неизвестен, использовать формулу 
(2) затруднительно. Построим для и„ функцию ра- 
спределения, предположив, что она относится к 
классу $Т$-распределений с некоторыми параме- 
трами. 

Определение. Пусть 

ё в О) = ^ } ехр [іх(р -0] х 


= <Р(Ф~'(Р,))іёв( а )У\ 

а 2 =ср{ф-\р 2 ))-[ ёѳ т\ 

«і = а-а і Ф\р 1 ), а 2 = Ъ + а 2 Ф \р 2 ), (5) 

где ср(х ), Ф(х) - плотность и функция распределе- 
ния стандартной нормальной случайной величины 
соответственно. 

Доказательство теоремы основано на использо- 
вании равенств (3) и выполнении соотношений: 


хехр 




( шѴ 

\ 

сіі 

а 

р, = С к (іуіі = ф 

( \ 

а-а х 

Г2~Х 

/ 

—со 

к а 1 У 


- функция плотности а-устойчивого распределе- 
ния с вектором параметров (а, Д, с, р), где а - ха- 
рактеристическая экспонента, Д - коэффициент 
асимметрии, с - масштаб, р — среднее. Пусть /г,(х), 
/г 2 (х) - плотности нормального распределения с па- 
раметрами (й;,< 7; 2 ), /= 1,2. Пусть выбраны два дей- 
ствительных числа а, Ь, причем а<р<Ъ, и выполне- 
ны соотношения: 

\{а) = § ѳ {а), Л 2 (й) = я ѳ (й), 

а а оо со 

| /?,(х)г/х = | & ѳ (х)сІх, |й 2 (х)с/х = |§ в (хУх. (3) 

—со —со Ь Ь 

Назовем плотностью $Т$-распределения функ- 
циюДх) вида: 


Р 2 = \ К (?№ = 1 - Ф 


^ Ъ-а 2 ^ 


Для задания в (2) стандартной $Т$-распреде- 
ленной случайной величины в„ требуется знать 
первый и второй начальные моменты. Справедли- 
ва следующая теорема. 

Теорема 2. Пусть Е,~8(х,а,Ь,а,Р,с,р). Тогда пер- 
вый и второй начальный моменты вычисляются по 
формулам: 


Щ=аРх -оДФ "Ч аХРі +^(Ф“ 1 (д))] + 

ь 

+| Хё в (х) г/х + Ьр 2 + ст 2 [Ф ■ _1 (р 2 )р 2 + (р(ф ' (р 2 ))], 

а 


Дх) = 


/?, (х),х < а; 
Яѳ(х),хе[а,6]; 
й 2 (х),х >Ъ. 


(4) 


Математические ожидания и дисперсии (а,, а?), 
і= 1,2 однозначно определяются из равенств (3). 

х 

Поэтому 5(х) = | / (()<]( = 8(х,а,Ъ,а, Д,с, р) будет 


зависеть от шести параметров: а, Ъ , а, Д, с, р. Раз- 
личные возможные параметризации рассмотрены 
в монографии [13], где обсуждаются их преимуще- 
ства и недостатки. 

Обозначим через р х и р 2 вероятности отсечения: 


2 = + а, ) - с, (а + а, ) (р(Ф ( р х )) + 

ъ 

+|х 2 Я ѳ (х)г/х + (ст 2 + а\)р 2 +а 2 (а 2 +Ъ)(р(Ф~\р 2 )). 

а 


Доказательство теоремы основано на непосред- 
ственном интегрировании равенства (4). 

Проведем сравнение построенной по (3)-(5) 
плотности 8Т8 -распределения с плотностями нор- 
мального и а-устойчивого распределения с идентич- 
ными параметрами, для чего изобразим их на рис. 1 . 

Как следует из анализа рис. 1, плотности 
а-устойчивого и 8Т8 -распределений совпадают 
друг с другом при хе[а,Ь\. Кривая 2 отличается от 
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Рис.1. Плотности нормального, а-устойчивого и 5Т5-распределения при а=~0,1, Ь=0,1 (р^=рі=0,105), сс=1,5, р=0,1135, 
с=0,05, р=0, 00191: 1) плотность 5Т5-распределения; 2 ) плотность нормального распределения с математическим ожи- 
данием р=0, 00191 и дисперсией 0,004761; 3) плотность а-устойчивого распределения 



Рис. 2. Плотности 575-распределения при а=1,5, р=0,1135, с=0,05, р=0, 00191 и функция плотности нормального распределе- 
ния с математическим ожиданием р и дисперсией 0,00476: 1) плотность нормального вероятностного закона; 2) плот- 
ность 5Т5-распределения при а=~0,1, Ь=0,1, р ] =р 2 =0,105; 3) плотность 5Т5 при а=~0,2, Ь=0,2, р ] =р 2 =0,031; 4) плот- 
ность 5Т5-распределения при а=~0,3, Ь=0,3, р\=р г =0,015 
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кривых 1, 3 вследствие ненулевого коэффициента 
асимметрии /3. Кроме того, характеристическая эк- 
спонента а не равна двум (случай нормального ве- 
роятностного закона). Далее, для х>тах{|а|,|й|} хво- 
сты кривой 1 находятся между кривыми 2 и 3. Это 
говорит о том, что редкие события при их описа- 
нии с помощью 8Т8 -распределения происходят с 
более высокой вероятностью, чем при использова- 
нии нормального закона. 

Зависимость плотности 8Т 3 - рас предел е н и я от 
параметров а, Ь приведена на рис. 2. Как следует из 
анализа рис. 2, с ростом абсолютного значения 
уровня отсечения а толщина хвостов падает, а 
функция плотности концентрирует свои значения 
около моды /и. Кроме того, при определенных зна- 
чениях параметров плотность 8Т8-распределения 
может иметь более толстые хвосты, чем плотности 
нормального и «-устойчивого распределений. По- 
этому можно утверждать, что построенная в соот- 
ветствии с (3)-(5) функция Дх) обладает высокой 
адаптивной способностью к описанию эмпириче- 
ских данных, достигаемой варьированием шести 
параметров. 


3. Эконометрический анализ данных 


Используем построенное выше 5Т$- рас преде- 
лен не для модификации ОАКСН(1,1). Предполо- 
жим, что в выражении (1) и~8{х,а,Ъ,а,Р,с,іл). Вы- 
берем стандартную 8Т8-распределенную случай- 
ную величину г„~5(х;-5,92;3,33;1,85;0,6;-0,1;0), 
имеющую нулевое математическое ожидание и 
единичную дисперсию. 

Предполагая независимость логарифмических 
приращений, применим для оценивания коэффи- 
циентов у,аф модели (1) метод максимального 
правдоподобия (вопросы применимости метода в 
случае условных плотностей подробно рассмотре- 
ны в [14]) и вычислим максимум функции Ь: 

т 

/=1 

или функции 1п Ь: 

т 

іпі=2/„ (6) 

і=і 

где /,'=Дх,|/^_і) - функции условной плотности 8Т8- 
распределения, определенные равенством (4), т - 
число наблюдений, I - функция правдоподобия. 

Поиск максимума (6) осуществляется в соответ- 
ствии с выполнением необходимого условия суще- 
ствования экстремума функции трех переменных: 


дЬ 

дсо 


= 0, 


дЬ 

да 


= 0 , ^4 = 0 . 

др 


( 7 ) 


Решение нелинейной системы (7) в предполо- 
жении единственности экстремума в некоторой 
расчетной области может проводиться любым ите- 
рационным методом: методом наискорейшего спу- 
ска, сопряженных градиентов и т. п. 


После оценивания коэффициентов а>,а,р и под- 
становки их в (1) остается провести статистическое ис- 
следование надежности предложенного метода 
ОАКСЩІД) при прогнозировании волатильности. 
Для этого используем хорошо известную статистику 
Льюнга-Бокса проверки гипотезы Н 0 о равенстве нулю 
первых т автокорреляций [15], где т<п. Так как при 
условии существования четвертого начального момен- 
та Т|й„|<со процесс ОАКСЩІД) может быть записан в 
виде АКМА процесса с параметрами р= 1 и д= 1 , то есте- 
ственно рассмотреть нормированную выборочную ав- 
токорреляционную функцию остатков а„ вида 

л П ~ к /V л / я л 2 

Гк = ^сцсц+к /2^ я/ ,к = 1,2,3,... 

і=і / /=і 

Далее, запишем статистику 

л2 

т Гк 

У=п(п + 2)]Г 

которая, как известно, при достаточно больших п 
будет иметь ^-распределение с {т-р-ц) степеня- 
ми свободы, если теоретические значения параме- 
тров модели (1) неизвестны. Наконец, вычисление 
г к следует проводить для остатков й„=и „ 2 и 
а„=(и п —а)/(тІ соответственно до и после примене- 
ния ОАКСЩІД). 

Гипотеза # 0 принимается, если у<Хл-к т ~Р~4)ь г Д е 
5 - уровень значимости критерия, и отвергается в про- 
тивном случае. Соответственно, ОАКСН(ІД), опре- 
деляемый выражением (1) с коэффициентами, удо- 
влетворяющими (7), является статистически надеж- 
ным с уровнем значимости 5, если у<Хл-к т ~Р~4)- 

Для моделирования вероятных в будущем зна- 
чений временного ряда Н„ вместе с вычислением 
статистики Льюнга-Бокса потребуется проверить, 
будут ли приращения и„, найденные в соответствии 
с (2), иметь 8Т8-распределение. Подгонка эмпири- 
ческого распределения к 8Т8 по выборочным дан- 
ным может осуществляться различными методами, 
например, подгонкой характеристической функ- 
ции или аппроксимацией плотности с использова- 
нием быстрого преобразования Фурье. В данной 
работе использовалась запись плотности 8Т8-ра- 
спределения через интеграл Золотарева [13] с его 
последующим численным интегрированием ква- 
дратурной формулой Симпсона. Качество подгон- 
ки проверялось с помощью ^ 2 -теста [14]. 

4. Анализ эмпирических данных 

Применим метод 8Т8-ОАКСЩ1Д) для модели- 
рования временного ряда к„ долларовых котировок 
акций РАО ЕЭС РТС, для чего используем 371 значе- 
ние за период с 04 января 2003 г. по 30 июня 2004 г. 

Отметим, что і/ н =1п(/? й ) — 1п(/г„_,), л=1,2,..., обла- 
дают следующими параметрами: среднее - 1 ,9- 10 3 , 
дисперсия - 8,2- 10 4 , третий центральный момент 
- 2,6- 10 6 , четвертый - 4,1 -ІО' 6 , коэффициент 
асимметрии - 1,1- 10 ‘, куртосис (эксцесс) - 6,15. 
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Проведенные численные расчеты показали, что 
и~8(х,а,Ь,а,р,с,/и ) с вероятностью 0,99, причем 
й=-4-10 -2 , *=4,5-10 2 , а=1,9; >8=1,110 ‘, с=1,8-10 2 , 
и= 1 ,9- 10“ 3 . При этом значение / 2 -теста было равно 
21,67. 

Дневная волатильность <т„ в (1) вычислялась по 
последним к наблюдениям: 

= т^тХ(Ч,.-й) 2 , 

к А /=1 


где и = 


1 і 

— Х м и-< _ выборочное среднее, к - лаг (за- 

к ы 


держка) временного ряда. Коэффициенты модели 
(1) были оценены методом максимального правдо- 
подобия с логарифмической функцией правдопо- 
добия (6). Система (7) разрешалась методом на- 
искорейшего спуска с погрешностью е=10 3 . На- 
чальное приближение выбиралось нулевым. 

Величины параметров $Т8-ОАКСН(1,1) при 
различных лагах к приведены в таблице. 


Таблица. Значения коэффициентов 5Т5-САРСН (1,1) при 
различных лагах 


к 

СО 

а 

Р 

7 

V 

2 

ю- 4 

0,20 

0,77 

0,03 

з-ю- 

5 

ю- 

0,21 

0,78 

ю- 

ю- 

10 

ю- 

0,16 

0,83 

ю- 

ю- 


Приведем значения статистики Льюнга-Бокса 
у до и после применения САКСН(ІД). Так, при к = 3 
до и после применения она была равна 29,58 и 0,09 
соответственно, при к= 5 - 48,6 и 1,11, при &=10 - 
60,13 и 6,78. Пороговые значения распределения 
х\-і т ~ 2) при уровне значимости 5=0,05 и т= 3, 
т= 5, т= 10 равны 3,8414; 7,8147; 15,5073. Таким об- 
разом, метод 8Т8-САКСН(1,1), построенный на ос- 
нове (1) с коэффициентами из таблицы, является 
статистически надежным с вероятностью 0,95. 


Применим §Т$-ОАКСН(1,1) для имитационно- 
го моделирования долларовых котировок акций 
РАО ЕЭС РТС. При этом используем найденные ра- 
нее волатильности сг„ и сгенерируем случайную по- 
следовательность е,={у„}„ 3 ™ 0 , имеющую 8Т8-распре- 
деление с нулевым средним и дисперсией единица: 

е п ~5(х;~ 5,92; 3,33;1,85; 0,6; -0,1; 0). 

Задание значений в„ проводилось согласно 
классической схеме [17]: на интервале [0,1] генери- 
ровалась последовательность равномерно распре- 
деленных случайных величин после чего при 
фиксированном п решалось трансцендентное ура- 
внение относительно у„ : 

х п=$(уЛ ( 8 ) 

где 5(у) = 5(у; - 5, 92;3, 33; 1, 85; 0, 6; -0, 1; 0) = Г /( I) сП 

—00 

- функция стандартного 8Т8-распределения. Ре- 
шение уравнения (8) осуществлялось методом ка- 
сательных с нулевым начальным приближением и 
точностью 10“ 4 . Относительная погрешность вычи- 
сления корней (8) не превосходила КТ 6 . 

Последовательность е„ была использована в (2) 
для определения логарифмов приращений и„ и для 
вычисления вероятных значений долларовых 

котировок акций РАО ЕЭС РТС: 

•ехр(ст„ +1 г„ +1 + а„ +1 ). 

Полученные /г„ т , рассчитанные с лагом 5, и ис- 
ходные исторические данные к п приведены на 
рис. 3. Кроме того, на рис. 3 изображена относи- 
тельная погрешность <5„ между ними: 


Как следует из анализа рис. 3, относительная 
погрешность 5„ не превосходила 5 %. Максимум 8„ 
достигался на 206 день торгов и был равен 4,8 %. 
Средняя величина погрешности составила 2,28 %. 



Рис. 3. Сравнение исторических и модельных значений котировок акций РАО ЕЭС РТС за период с 04 января 2003 г. по 30 ию- 
ня 2004 г.: 1) временной ряд исторических данных, 2) временной ряд, рассчитанный 5Т5-ОАРСН(1,1) с лагом 5, 3) отно- 
сительная погрешность, % 
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Согласно рис. 3, наблюдаются три всплеска 
значений <5„: с 29.04.03 по 30.05.03, 20.10.03 и с 
14.04.04 по 17.05.04. Это связано с резкими скачка- 
ми котировок акций в данные промежутки време- 
ни, так как волатильность цен составляла от 7 до 
15 % в день. Отметим, что такие резкие изменения 
цен являются одной из причин отказа от допуще- 
ния о нормальном распределении е„ при вычисле- 
нии и„ в выражении (2). Как показывают расчеты, 
если вместо г„~5'(х;-5,92;3,33;1,85;0,6;-0,1;0) взять 
е„~/Ѵ(0,1), то максимум относительной погрешно- 
сти будет равен 17,1 % при средней величине по- 
грешности 7,3 %. Вышесказанное дает основание 
сделать вывод о том, что 8Т8-ОАКСН(1,1) удовле- 
творительно описывает исходные данные и позво- 
ляет моделировать их значения с небольшой по- 
грешностью. 
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5. Заключение 

На основе одномерного 8Т8 -распределения по- 
строен модифицированный эконометрический ме- 
тод 8Т8-ОАКСН(1,1). Найдены аналитические вы- 
ражения для первого и второго начальных моментов 
8Т8-распределения. Неизвестные параметры моде- 
ли 8Т8-ОАКСН( 1,1) определены численно методом 
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Льюнга-Бокса осуществлено статистическое иссле- 
дование надежности предложенного эконометриче- 
ского алгоритма. Далее, метод 8Т8-ОАКСН(1,1) 
был применен для имитационного моделирования 
цен обыкновенных акций РАО ЕЭС РТС с лагом 5. 
Показана более высокая точность вычислений по 
сравнению с классической моделью САКСН(і,І), в 
которой делается допущение о нормальности лога- 
рифмов дневных приращений. 
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